PROCEEDINGS 

OF SCIENCE 



Nucleon and /V ( 1535) Distribution Amplitudes 



V. M. Braun a , S. Collins M. Gockeler C. Hagen a , R. Horsley b , Y. Nakamura 

D. Pleiter d , P. E. L. Rakow e , A. Schafer R. W. Schiel* G. Schierholz /, H. Stiiben 

s, J. M. Zanotti b 

" Institutfiir Theoretische Physik, Universitat Regensburg, 93040 Regensburg, Germany 

b School of Physics and Astronomy, University of Edinburgh, Edinburgh EH9 3JZ, UK 

c Center for Computational Sciences, University ofTsukuba, Tsukuba, Ibaraki 305-8577, Japan 

John von Neumann Institute NIC / DESY Zeuthen, 15738 Zeuthen, Germany 
e Theoretical Physics Division, Department of Mathematical Sciences, University of Liverpool, 
Liverpool L69 3BX, UK 
f Deutsches Elektronen-Synchrotron DESY, 22603 Hamburg, Germany 
8 Konrad-Zuse-Zentrum filr Informationstechnik Berlin, 14195 Berlin, Germany 
E-mail: 



rainer . schiel@physik . uni-regensburg . de 



(The QCDSF Collaboration) 



The QCDSF collaboration has investigated the distribution amplitudes and wavefunction nor- 
malization constants of the nucleon and its parity partner, the A^*(1535). We report on recent 
progress in the calculation of these quantities on configurations with two dynamical flavors of 
&{a) -improved Wilson fermions. New data at pion masses of approximately 270 MeV helps in 
significantly reducing errors in the extrapolation to the physical point. 



The XXVIII International Symposium on Lattice Field Theory 
June 14-19,2010 
Villasimius, Sardinia Italy 

* Speaker. 



© Copyright owned by the author(s) under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike Licence. 



http: //pos . sissa. it/ 



Nucleon and N* (1535) Distribution Amplitudes 



R. W. Schiel 



1. Introduction 

Distribution amplitudes (DAs) provide access to hadron wave functions and can be related to 
form factor data through perturbative QCD and light cone sum rules. This makes DAs interesting 
quantities to be investigated within the framework of lattice QCD. Distribution amplitudes for the 
nucleon and its parity-partner, the TV* (1535), have been calculated by the QCDSF collaboration in 
[jl|, ||] on lattices with Nf = 2 degenerate Clover fermions and pion masses of 400 to 1400 MeV. In 
this update, we present new data obtained for a 32 3 x 64-lattice with j8 = 5.29, K = 0.13632 and 
270 MeV pions. This new data is significantly closer to the physical point than before and will help 
us in reducing the uncertainties in the extrapolation to the physical point. 

The nucleon distribution amplitudes are defined as follows [Q]. In the infinite momentum 
frame, with transverse momentum components integrated out and only leading twist components 
considered, the nucleon wave function for the three-quark Fock state can be written as 

\N,t) =f N I ^fkL{\uHx l )uHx2)dHx 3 ))-\uHx l )dHx 2 )uHx 3 ))} (1.1) 
J 2y24xiX2X3 

where xt are the longitudinal momentum fractions, the arrows indicate the nucleon and quark spins, 
f[dx] = Jq 1 dx\dx2dx 3 b{ \ —x\—X2 — x 3 ), fy is the leading-twist normalization constant and <p(xj) 
is the nucleon distribution amplitude. We also calculate the normalization constants for the next- 
to-leading twist wave functions, X\ and X 2 \ see Refs. [0, § for the definition. 



2. Lattice procedure 

On the lattice, only moments of the distribution amplitude, 

(p lmn = y n ^ x ] x ' lX »VJ<p(x 1 ,X2,X 3 ) 

can be accessed. And even this is limited to the first and second moments, since for higher mo- 
ments, mixing with lower dimensional operators occurs, which would make the calculation in- 
creasingly difficult. Moreover, for higher moments, the matrix elements have to be evaluated at 
higher momenta which leads to noisier data. As can be seen in the following, even for the second 
moments we have barely enough statistics to get reasonably small error bars, and third moments 
would be much worse than that. 

It is useful to expand the wave function in multiplicatively renormalizable terms [Q] : 

? f — — 

(pfain ) = 120x1x2x3-^ 1+cio(jci -2x 2 + x 3 )L 3 ' s o +cu(xi -x 3 )L 9 ft> 

+c 20 [1 + 7(jc 2 - 2x { x 3 - 2x 2 2 )] L 3Pa +c 2l (1 -4x 2 ) (jcj -x 3 )L 9 A> 
+c 22 [3 - 9x 2 + %x\ - 12*1*3] L$ + ■ ■ •} 

where L = a s (n) / a s (Ho) and c, ; - are the so-called "shape parameters". The c ; y are given by lin- 
ear combinations of the moments of the distribution amplitude (p lnm with I + m + n < i. Unlike 
the moments, which are constrained by the momentum conservation condition (p ,mn = <p(' +1 ) m " + 
yltjn+vjn _|_ (plmin+i)^ ^ e we independent nonperturbative parameters. In practice, they are ob- 
tained from the data on (p Imn using constrained fits. 
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AT* (1535) (right, green line). 
Figure 1: Effective mass plot (a) and dispersion relation (b). Only statistical errors are shown. 



The moments of the distribution amplitudes are calculated through matrix elements of the form 

{@{x) aM J ? {y) x ) 

where jV is a smeared nucleon interpolator and is a local three-quark operator with up to two 
derivatives. The operators & that are used to calculate the DAs have been classified according to 
irreducible representations of the lattice symmetry group in order to avoid mixing of operators [j|] 
and they have been non-perturbatively renormalized [Hj. 

Separation of the parity -plus (nucleon) and parity-minus (TV* (1535)) states has been achieved 
by using the generalized Lee-Leinweber parity projector (the non-generalized version of which 
was introduced in [0]), ^(1 + ^74), where Y is a suitable product of y,- matrices and m and E are 
the mass and energy, respectively, of the baryon. Please note that this projector might have to be 
modified for certain operators and momenta, but it is perfectly suited for the operator-momentum 
combinations that are used here. In Fig. [h] one can see that a reasonable mass splitting is achieved 
between the nucleon and TV* (1535) at momentum p 2 = 0. At higher momenta, the data becomes 
noisier, but it is still possible to find plateaus for both the nucleon and TV* (1535). As a consistency 
check, the dispersion relation has been looked at and it is fulfilled within error bars (see Fig. |lb| ). 

3. Results 

3.1 Normalization Constants 

The results for the normalization constants fy, /n* and At 2 are shown in Fig. ^[ In comparison 
to the plots of the nucleon normalization constants shown in [|l|], adding the new data point at 270 
MeV pion mass makes a large difference, reducing the distance for the chiral extrapolation by 
more than a factor of two. Thus we can avoid using data with pion masses > 1 GeV which strongly 
influenced our previous results [p. 
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Figure 2: Dependence of fy and fy* (left) and Ai and A2 (right) on the pion mass. Only statistical 
errors are shown. 

It is interesting to note that on the lattice with 270 MeV pion mass, it is the first time to 
see a significant difference between the leading-twist normalization constants of the nucleon and 
AT* (1535). While fy and fy* agree - within error bars - on the lattices with heavier pions, this is 
no longer true for the new data point. However, we do not yet have an explanation why this is the 
case or whether it should have been expected. 

3.2 Shape parameters 

The shape parameters are calculated from the moments of the DAs as described above. Since 
the shape parameters describe the deviation from a symmetric momentum distribution, their relativ 
error is larger than the error of the moments of the DAs. Also, as can clearly be seen in Fig. ||, 
the second order shape parameters C2j are noisier than the first order ones, cy. This can be mostly 
attributed to the fact that the calculation of the C2j involves matrix elements at p 2 = 2, while one 
needs only the cleaner p 2 = 1 matrix elements to calculate the c\j. Also, the nucleon error bars are 
smaller than the 7V*(1535)'s, which is not surprising given that the AT* (1535) has a larger mass and 
thus a shorter plateau than the nucleon. 

The first order shape parameters show a significant difference between the nucleon and its 
parity partner. This difference is also visualized in the barycentric plot, Fig. f|: The AT* (1535) 
wavefunction is much more strongly peaked towards the spin-carrying quark (quark label "1" in 



Eq. 1.1) than the nucleon. 

The second order shape parameters still have very large error bars and are thus not included 
in the barycentric plot. However, we hope that increased statistics will reduce these error bars to a 
reasonable size. While C20 is compatible with both for the nucleon and its parity partner, C21 and 
C22 tend to be negative in both cases. 

4. Conclusions and Outlook 

Nucleon and AT*( 1535) distribution amplitudes and the wavefunction normalization constants 
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Figure 3: Shape parameters for the nucleon and N* (1535) on the 32 3 x 64 lattice at 270 MeV pion 
mass. Only statistical errors are shown. While the error bars for the first moments are fine, more 
statistics will improve the errors bars for the second moments. 




Figure 4: Barycentric plot of the nucleon (left) and Af*(1535) (right) distribution amplitudes, cal- 
culated for the central values of the shape parameters. Due to the large error bars on the second 
moments, only the first moments have been used for this plot. 



/iv> In* and X\2 nave been calculated at pion masses as low as 270 MeV. Unlike for higher pion 
mass lattices, a difference between the nucleon and TV* (1535) wavefunction normalization con- 
stants can now be seen. The distribution amplitudes for the nucleon and its parity partner are also 
significantly different, with the 7Y*(1535)'s wavefunction being more peaked towards the spin- 
carrying quark than the nucleon's wavefunction. 

To improve the accuracy of the results on the nucleon and AT* (1535) DAs, more statistics are 
being collected for the 32 3 x 64 lattice with 270 MeV pions. Also, calculations at the same pion 
mass but a larger lattice volume will be performed to get a handle on finite volume effects. Finally, 
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the data will have to be extrapolated to the physical point. 
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